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SECTION  1 


INTRODUCTION 


This  report  describes  the  present  status  of  a  computer  program  for  sim¬ 
ulating  the  RF  behavior  of  space-based  radar  (SBR)  phased  array  lens  antennas . 
A  typical  SBR  candidate  lens  is  shown  in  Figure  1.  The  transmitter  provides 
the  space  feed  that  illuminates  the  lens  during  radar  transmit.  The  space 
feed  may  be  in  the  form  of  several  independently  controlled  beams  that  pro¬ 
vide  lens  radiation  pattern  shaping,  adaptive  nulling,  and  time  delayed  lens 
sector  illumination.  (The  time  delayed  compensation  is  required  because  of 
significantly  different  transmitter-to-lens  path  lengths  to  different  points 
on  the  lens).  Beam  steering  and,  possibly,  amplification  are  performed  by 
the  lens. 

A  lens  comprising  two  arrays  of  thin  wire  radiators  separated  by  a 
ground  screen  currently  is  modeled  in  the  simulator.  The  theory,  with 
straightforward  modifications,  can  be  extended  to  apply  to  other  lens  types 
such  as  those  composed  of  microstrip  arrays. 

A  typical  wire  radiator  lens  section,  that  can  be  analyzed  by  the  sim¬ 
ulator,  is  shown  in  Figure  2.  For  simplicity,  dipole  radiating  elements  are 
indicated.  The  simulator  also  is  capable  of  analyzing  other  element  types; 
e.g.,  folded  dipole,  turnstile,  and  parasitic.  The  array  lattices  may  be 
skewed.  Focusing  and  scanning  of  the  lens-transmitted  main  beam  is  accom¬ 
plished  by  electronic  modules  interconnecting  the  radiating  elements  between 
the  illuminated  and  nonilluminated  arrays  (Figure  3).  These  modules  also 
may  provide  power  amplification. 

A  successive  approximation  method  is  employed  whereby  infinite  array 
analysis  provides  a  first  approximation  of  array  currents  and  patterns. 

Higher  approximations  then  are  obtained,  with  relative  ease,  by  methods  that 
draw  upon  these  first  approximation  currents.  The  higher  approximations  pro¬ 
vide  corrections  to  the  currents  on  elements  in  the  vicinity  of  discontinui¬ 
ties  in  periodicities;  e.g.,  lens  edges,  lens  section  interfaces,  and  failed 
modules.  A  section  of  uniform  periodicity  of  a  lens  is  called  a  cell. 

The  first  approximation  patterns,  although  involving  large  numbers  of 
elements,  are  computed  from  closed  form  expressions.  Since  the  number  of 
elements  with  currents  perturbed  from  their  infinite  array  values  is  expected 
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to  be  relatively  small,  higher  approximation  "correction  patterns"  may  be 
computed  by  direct  summation  of  the  field  from  each  element.  The  superposi¬ 
tion  of  the  first  approximation  and  correction  patterns  results  in  the  over¬ 
all  pattern.  This  approach  takes  maximum  advantage  of  the  nominal  periodici¬ 
ties  and  largeness  of  the  arrays.  It  is  as  computationally  manageable  in 
the  radar  receive  mode  as  in  the  radar  transmit  mode  without  resorting  to 
reciprocity.  This  is  important  since  the  modules  may  be  nonreciprocal.  This 
approach  also  may  avoid  the  need  for  a  Fast  Fourier  Transform  (FFT)  for  pat¬ 
tern  computation.  An  FFT  is  not  practical  if,  for  example,  fine  angular 
pattern  resolution  in  the  vicinity  of  a  null  is  required. 

An  infinite  array  analysis  usually  assumes  plane  wave  excitation.  Dur¬ 
ing  radar  transmit,  the  illuminating  field  is  expected  to  differ  smoothly 
from  a  plane  wave  across  the  face  of  the  lens.  This  difference  is  taken  into 
account  by  expanding  the  beamformer  field  illuminating  a  cell  in  plan  waves 
and  analyzing  the  cell  separately  for  each  plane  wave  component. 

The  first  approximation  method  is  an  extension  of  a  moment  method, 
plane  wave  expansion  technique  recently  suggested  by  Munk  and  Burrell  [1]. 

It  is  applicable  to  lens  arrays  composed  of  arbitrarily  shaped  wire  radiating 
elements.  The  elements  may  be  inclined  with  respect  to  the  array  plane. 
Feedline  scattering,  a  nonflat  lens,  amplitude  and  phase  adjust  modules  be¬ 
tween  the  arrays,  and  an  imperfect  ground  screen  also  are  taken  into  account. 
(The  nonflat  lens  is  modeled  in  "piecewise"  fashion  by  allowing  the  cells  to 
be  tilted  with  respect  to  each  other.)  The  plane  wave  expansion  technique 
facilitates  computing  the  array-to-array  coupling  through  the  imperfect 
ground  screen  and  suggests  straightforward  extensions,  not  discussed  here, 
applicable  to  dielectric  support  sheets  (e.g.,  Kevlar  and  Kapton)  on  which 
the  radiators  may  be  mounted. 

The  infinite  array  lens  analysis  facilitates  computation  by  providing 
the  first  approximation  solution  to  a  finite  array  lens.  It  also  provides 
clues  to  the  performance  of  several  important  finite  (but  large)  lens  char¬ 
acteristics.  For  example,  the  array  active  impedance  variation  with  module 
phase  setting,  an  important  parameter  for  module  designers,  can  be  assessed 
directly.  The  effect  of  imperfect  ground  screen  "feedback"  from  the  target 
side  to  the  feed  side  of  the  lens  also  can  be  observed,  and  lens  warping 
(nonflat  lens)  treated. 
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The  application  and  generation  of  infinite  array  theory  and  the  moment 
method  in  arriving  at  a  first  approximation  to  the  lens  radiating  element 
currents  has  been  detailed  in  [I]-  This  theory  has  since  been  improved  and 
corrected  and  the  development  consolidated.  The  revision  is  given  in 
Section  ?.  The  principal  improvement  lies  in  replacing  a  pulse  expansion- 
pulse  matching  moment  method  with  a  pulse  expans  ion  ■  y  'i:.t  aa: ching-f inite  dif¬ 
ference  operator  moment  method.  Section  4  shows,  through  comparison  with 
another  method,  that  i he  "finite  difference"  method  is  computationally  more 
accurate  and  efficient  than  the  "pulse  matching’'  method.  .Also,  radiator 
arms  inclined  with  respect  to  the  array  plane  are  modeled  far  simply  with 
the  "finite  difference"  method. 

The  theory  underlying  the  higher  approximations  current  computation 
and  the  first  and  bight r  approximations  pattern  computation  is  essentially 
unchanged  from  that  in  [5]  and  not  repeated,  hero. 

The  main  computer  programs  and  subprograms  are  described  in  Section  3. 

A  general  module  model  in  terms  of  two-port  scattering  parameters  is  not  yet 
implemented  but  should  be  available  shortly.  Two  limited  module  models  -  an 
S  parameter  transmission  line  passive  phase  shifter  model  and  the  impedance 
representation  described  in  [5]  -  presently  are  implemented  and  are  included 
in  the  program  descriptions.  Data  input  is  described  in  Section  3.3.11. 

Section  4  gives  results  of  applying  the  first  approximation  plane  wave 
expansion  (infinite  array)  moment  method  to  active  impedance  and  current  dis¬ 
tribution  computations  for  dipole  array  elements  with  feedline  scattering 
accounted  for.  Results  for  "swept  baek'-dipole  arrays  also  are  given.  This 
is  the  first  time  theoretical  act Lve  impedance  and  current  computations  of 
infinite  planar  arrays  of  radiating  elements  inc1  ined  witn  respect  to  the 
array  plane  have  been  reported. 


SECTION  2 


IMPROVEMENT  AND  CONSOLIDATION  OF  FIRST  APPROXIMATION  THEORY 


A  combined  plane  wave  expansion  moment  method  technique  is  applied  in  ob¬ 
taining  a  first  approximation  of  the  radiating  element  currents  in  a  lens.  It  is 
applicable  to  lens  sections  (cells)  composed  of  many  identical  arbitrarily 
bent  wire  radiating  elements  in  each  of  two  parallel  arrays.  It  also  accounts 
for  amplitude  and  phase  adjust  modules  interconnecting  the  arrays  and  an  im¬ 
perfect  ground  screen  between  the  arrays.  In  particular,  the  technique  per¬ 
mits  the  modules  to  be  progressively  phased  although  the  module  amplitude  ad¬ 
justments  must  be  uniform  throughout  each  cell. 

A  lens  formed  by  two  planar,  periodic  infinite  dipole  arrays,  a_  and  on 
either  side  of  a  finitely  conducting  infinite  ground  screen  is  shown  in  Fig- 

->0X 

ure  2 .  An  exciting  plane  wave,  E  ,  that  is  arriving  from  the  negative  z  side  is 
directly  incident  on  array  a  (Figure  3).  Amplitude  and  phase  adjust  modules 
connect  array  a  to  array  b^.  Array  a  faces  the  radar  transmitter  during  radar 
transmit  and  faces  the  target  during  receive.  Array  a_  thus  is  always  on  the 
illuminated  side  and  array  b^  is  always  on  the  nonilluminated  side.  Each  array 
lattice  may  be  general  rectilinear  (skewed)  of  which  "rectangular"  is  a  special 
case.  Array  a  elements  may  differ  from  array  b^  elements;  however,  each  element 
within  an  array  is  assumed  identical. 

The  following  analysis  of  this  lens  is  not  restricted  to  dipole  array 
elements.  Each  element  can  be  a  collection  of  thin  bent  wires;  e.g.,  folded 
dipole  or  dipole  with  parasitic  scatterers.  The  elements  can  be  inclined  with 
respect  to  the  array  plane.  Also  feedline  scattering  is  considered. 

The  analysis  is  structured  around  a  port  representation.  Each  pair  of 
terminals  entering  an  array  element  is  assumed  to  be  a  port.  The  port  vol¬ 
tages  and  currents  are  shown  in  Figure  3.  The  superscript  (a  or  b)  indicates 
the  array.  The  subscript  indicates  the  element  location  with  reference  to 
Figure  4  where  mn  denotes  the  m^  element  along  a  line  parallel  to  the  lattice 
e  coordinate  and  n1"*1  element  along  a  line  parallel  to  the  lattice  y  coordinate. 
The  mn  »  00  elements  lie  on  the  z  axis  and  are  referred  to  as  the  reference 
elements. 
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Balanced  mode  transmission  line  theory,  applied  to  the  mn  elements' 
feedlines  and  module,  is  employed  in  relating  the  mntb  element  array  a_  port 
current  and  voltage  to  the  mnt'ri  element  array  port  current  and  voltage. 

This  relationship  is  combined  with  a  second  relationship  between  the  port  vol¬ 
tages  and  currents  in  achieving  a  solution  for  all  port  voltages  and  radiating 
element  current  distributions.  The  second  relationship  employs  a  plane  wave 
expansion  moment  method  technique  that  accounts  for  all  radiation  coupling 
including  field  penetration  through  the  ground  screen  and  feedline  (unbalanced 
mode)  scattering.  The  feedlines,  whether  coaxial  or  twin  wire,  are  modeled 
(Figure  5)  as  single  thin  wires.  If  the  actual  feedline  is  twin  wire,  the 
radius  of  the  feedline  model  wire  is  chosen  such  that  the  model  and  the  actual 
twin  wire  exhibit  approximately  the  same  scattering. 

The  Figure  5  model  assumes  zero  unbalanced  mode  feedline  currents  at  the 
element  ports.  The  model  also  does  not  allow  unbalanced  feedline  currents  to 
flow  through  the  modules.  These  assumptions  can  be  removed  by  employing  more 
general,  but  more  complicated,  models,  e.g.,  treating  modules  as  four-ports 
instead  of  two-ports  and  permitting  multiple  wire  junctions  at  the  feed  ports. 

2.1  PORT  REPRESENTATION  AND  SOLUTION 

Array  a  and  array  b^  port  currents  are  determined  by  requiring  that  they 
satisify  the  port  boundary  conditions* 


mn 


vaa  +  vab  +  vex 

ran  ran  ran 


a) 


mn  mn 


(2) 


where 


&  b  3  b 

V  and  V  ,  the  total  port  voltages,  are  related  to  I  and  I  via 

ran  mn  ran  ran 

mn^  module  (including  feedlines)  two-port  parameters; 

aa  th  i^ex 

V  is  the  array  a  mn  element  port  voltage  with  E  =0,  the  ground 

ran  — 

screen  assumed  perfectly  conducting,  and  the  array  a_  ports 
excited  with  the  I  as  ideal  current  sources  (Figure  6); 


uhe 


*  A  suppressed  expfjwt]  temporal  variation  is  assumed  throughout  the 

development. 
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Y 


PERFECT  CONDUCTOR 


Figure  6.  Problem  Defining  V 


aa 

mn 


SSSS\__  ^ -  PERFECT  CONDUCTOR 

PX 

Figure  3.  Problem  Definino  V _ 


bb  .  aa 

V  is  the  array  b  counterpart  to  V  ; 

mn  —  mn 

Vab  is  the  array  a  mnth  element  port  voltage  with  Eex  =  0,  all 
mn  — 

b 

array  a  ports  open-circuited ,  and  the  I  as  ideal  current 
—  mn 

sources  exciting  the  array  b  ports  (Figure  7); 

6X  til 

V  is  the  array  a  mn  element  port  voltage  with  all  array  a 

ports  open-circuited  and  the  ground  screen  assumed  perfectly 

conducting  (Figure  8). 

In  arriving  at  (1)  and  (2),  the  finite  conductivity  of  the  ground  screen 
is  considered  significant  with  regard  to  feedback  from  array  b  to  array  a . 
Severe  pattern  degradation  may  result  if  the  module  amplification  is  compara¬ 
ble  to  the  ground  screen  attenuation;  therefore,  it  is  important  to  preserve 
this  effect  in  the  model. 

Since  the  excitation  is  assumed  to  be  a  plane  wave,  the  relationship  of 
the  current  on  the  mnth  element  to  the  reference  (00th)  element  current  is 
known.  It  therefore  is  necessary  to  solve  (1)  and  (2)  for  the  reference  ele¬ 
ment  current  only.  Each  array  element  is  modeled  as  piecewise  linear  and  the 
array  is  viewed  as  a  collection  of  infinite  arrays,  each  one  corresponding  to 
one  linear  segment.  A  moment  method  procedure  is  applied  where  each  expansion 
function  is  one  of  the  infinite  arrays.  The  field  from  each  "expansion  array" 
is  expressed  as  a  doubly  infinite  sum  of  plane  waves.  One  advantage  of  this 
plane  wave  expansion  is  that  penetration  through  imperfect  ground  screens 
becomes  straightforward. 

The  method  presented  here  is  an  extension  of  Munk  and  Burrell's  work  in 
that  interconnecting  modules  are  introduced  between  arrays  and,  more  impor¬ 
tantly,  the  modules  may  differ  from  element  to  element.  In  particular,  the 

"v 

modules  may  introduce  a  linear  progressive  phase  and  thus  differ  in  a  wetl- 
defined  manner.  If  the  modules  are  identical,  the  incident  plane  wave 

EGX  (r)  =  E6*  (tf)  exp [ -jks  x  -jks  y  -jks  z  ] 

x  y  z 
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excites  port  currents  of  the  form 


=  100  eXP  l~  jkB,deSe  "  jkndySv  1 


1^  =  exp  [~  jkmd  s 

mn  00  r  e  e 


-jknd  s  ] 

y  y 


where 


— ►  A  A  A 

r  *  xx  +  yy  +  zz  is  a  field  point  position  vector; 


x,  y,  z  are  the  rectangular  coordinate  system  unit  vectors; 
k  is  the  propagation  constant; 

s  ,  s  ,  s  are  the  rectangular  coordinate  directional  cosines  of  the 
x  y  z 

propagation  direction  of  E**; 

dg  and  1^  are  the  lattice  e  and  y  interelement  spacings;  (Figure  4)  ; 

I*0  and  are  the  reference  element  port  currents. 

The  directional  cosine  s  is  given  in  terms  of  s  and  s  by  [  5  (Section  3.2.4)] 

e  x  y 

d  s  -Ays 
xx  y 


where  d^  and  Ay  are  defined  in  Figure  4. 

The  justification  for  (3)  and  (4)  can  be  demonstrated  from  linearity  con¬ 
siderations;  a  similar  form  applies  to  the  induced  port  voltages.  Equations 
(3)  and  (4)  are  not  sufficiently  general  if  the  modules  are  progressively 
phased.  If  the  phase  imparted  by  the  mn^  module  with  respect  to  the  refer¬ 
ence  module  is  -k(mde  ag  +  nd  a  ),  the  feedback  phenonmenon  due  to  the  imper¬ 
fect  ground  screen  and  impedance  mismatch  suggests  port  currents  of  the  forms 

I*  “  ^  I*.  .  exp  [-JJand  (s  -Km  )  _  jknd  (s  +ua  )  ]  (5) 

■»  u-_«  00  (u)  e  e  e  7  7  7  K  1 


1^  “  53  Inn/  \  exP  [-  Jkmd  (s  +uo  )  -  jknd  (b  +uo  )  ] 
u-_.  00(u)  e  e  e  7  7  7 
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where  ,  for  u  =  0, 

zaa  [s  ,  s  J  is  the  array  a  active  impedance  when  driven  by  port 
e  y 

currents  (  3  }  and  the  ground  screen  perfectly  conducting; 

zbb  [s  ,s  1  is  the  array  b  active  impedance  when  driven  by  port 
e  y 

currents  (  4  )  and  the  ground  screen  perfectly  conducting,- 

zab  [s  ,s  ]  is  the  active  mutual  impedance  from  array  b  ports  to 
e  y 

array  a  ports  with  array  b  excited  by  port  currents  (  4  )  and 

array  a  open-circuited  (feedlines  removed). 


U 


Specifically, 


taaTs  ,s  ] 
e  ’  yJ 


Vaa 

00 


00 


perfect  ground  screen 


bbr  t 

z  [se*  y1 


rbb 

00 

b 

00 


perfect  ground  screen 


,abtVsyl 


ab 

00 

b 

00 


imperfect  ground  screen,  array  a  ports 
opcn-circuitcd 


a  b 

where  I  and  XQ0  are  reference  element  port  currents  with  the  arrays 
excited  by  port  currents  related  by  {  3  )  and  (  4  ) ,  A  moment  method 

technique  for  computing  zaa[se»Sy],  zbb[se,s^],  and  zat>[se,sy]  is  des¬ 
cribed  in  Section  2.2.  Similar  relations  hold  for  u  /  0. 


A  means  for  determing  the  short-circuit  array  a  port  currents 

_ex  ex  r  ,  ...  i 

I  =  I  exp  f-  ]kmd  s  -jknd  s  J 
mn  00  1  eeJ  y  y 

,  -*ex 

due  to  E  with  the  ground  screen  perfectly  conducting,  is  described 

GX 

in  Section  2.2.  The  V  in  (  1  )  then  are  given  by  [5,  Eg. (3-20)] 

mn 


„ex  aar  ■,  ex  r  ,  , 

V  =  -  z  [s  ,s  J  I  „  exp  [-Tkrnd  s  -  -jknd  s  ] 
mn  1  e  yJ  00  *  1  J  e  e  J  y  yJ 


(7) 


ex  a 

where  the  reference  direction  for  I  is  the  same  for  I  in  Figure  5. 

mn  mn 

In  column  vector  notation,  (  5  )  and  (  6  )  can  be  represented 


as 


1 

mn 


X)  W)  exp  [-jkmd  (se  +  UV  '  jkndy(sy  +  UV  1 


(8) 
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where 


0 


s 

y 


]  i 


ex 

00 


Let  the  mn*'*1  module  (plus  feedlines)  scattering  parameter  matrix  be  denoted 

Si  b 

[S  ],  If  the  scattering  parameter  normalizing  impedances  are  R  and  R  for 
the  array  a_  and  sides  respectively, 


V 

urn 


[R]1/2([0]  +  [S  ])C 

ran  mn 


(10) 


i  = 
mn 


-{R]"1/2<(U]  -  f S  ] > C 

mn  mn 


{id 


where 


mn 
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a  b 

and  C  and  C  are  the  "incident"  scattering  variables.  Assume  that 
mn  mn 

C  may  be  expanded  in  "feedback”  modes,  as  is  j  ,  accordinq  to 
mn  mn 


5  =  7  c__/  x  exp  [- jkmd  (6  -in  a  )  -  jknd  (s  +ua  )  J 

mn  /  a  00  (u)  e'e  e  J  y  y  y ' 


Then  (  n  )  ,  (  10  )  ,  and  (  12  )  combine  to  vi.ei<! 


CD 

5ZeXP  t'j^e^e  'jkndyU“  1  UJ,/2(IU]  +  ls*”S00(u) 

u»-“  L 


^c^OOCu)  *  V00 


Another  equation  in  CQ0(u)  and  I  ,  obtained  from  (  >?  )  ,  (  11  ) 


and  (  12  ) ,  is 


^  ^  exp  L  -  j 


jkmd^ua^  -jknd^ua^  j  [R]  ([U]  -  [S^])^ 


mn  ;  00(u) 


+  ^OCu)  ~  0 


where  0  is  the  null  column  vector.  A  set  ol  simultaneous  equations  is  formed 

from  (13)  and  (.14).  The  modulo  "phase  set  tins-.;",  ■  and  u  ,  will  be  of  such 

t  y 

magnitude  that  [S  ’  repeats  every  M  x  j;  elements.  I'remult  iplication  of  (13) 
mn 

and  (14)  by  exp  [jl!"mvM'/M  +■  i?nvN'/"i,  where  v  is  an  integer,  summation  over 
one  "progressive  phase  period,"  l-rf-M  and  V-n'-N,  and  interchange  of  summations 
yields 
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N  M 

[F~'  ■  EE  exp  [  j2mn(v  -  u)M’/M  +  j27m(v  -  u)N’/N  ]  [ U] 

n*l  m»l  00) 


|MN[U]  u  =  v 
^  0  otherwise 


The  simplicat ions  resulting  in  the  right-hand  sides  of  (15)  and  (19)  assume 
that  (v-u)M'/M  and(v-u)N'/N  are  never  integers  simultaneously.  This  condition 
occurs  whenever  (v-u)  is  an  integer  multiple  of  M  and  N  simultaneously.  This 
condition  is  avoided  by  limiting  (15)  and  (19)  to  only  the  lower  valaed 
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equations  in  |vj  and  by  limiting  the  infinite  summations  to  only  lower  values 
of  |uj.  This  is  permissible  if  only  lower-ordered  ''feedback"  modes,  u,  are 
significant. 

A  special  case  of  interest  is  a  passive  lens  with  modules  replaced  by 
transmission  lines.  The  relative  lengths  of  line  determine  the  progressive 
phasing.  The  scattering  parameters  for  the  mn^  element  become 

!snJ  =  exp  t_j2nm  m"  'j2"n  7T  1  [W! 

if  R3  =  =  Rg,  where  is  the  characteristic  impedance  of  the  transmission 

line  and 


W  = 


expt-jkl 


exp[-jk£Q]  0 


where  ..  is  the  length  of  the  reference  element  transmission  line.  Equations 
(17)  and  (18)  reduce  to 


[F'  ]  = 
vu 


/T  MN  [U] 
/~R^  MN  (W  | 


V  -  u  -  0 
v  -  u  -  1  =  0 

otherwise 


(MN/kR  )  [U] 
f  I-’"  u  1  =  /  (-MN/ /r^  )  [W] 


v  -  u  *  0 

v  -  u  -  1  =  0 
otherwise 


Equations  (  15  )  and  (  16  )  then  become 


'  <'O0(v)  +  >R0  ^W'C00(v-l)  ^  7\ ^  1 0 0  ( v ) 


C  v  -  0 
0() 


v  =  U,i2,  . 
(20) 


-1-  ”■  -  — —  ( wlc  +  T  -  0  v  -  a .  f  1  * ,  .  .  .  (21 ) 

,R0  00  (v)  /^70  1  OO(v-l)  00  (v)  ' 
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Finally,  (  20  )  and  {  21  )  can  be  combined  to  form 


+  4  l*»1 


C  + 

C00(v) 


1 

Ro 


/R0  [u] 


] 


[W]c 


00(v-l) 


(  0  v  =  ±1, ±2,  .  .  .  (22) 

Equation  (  22  )  can  be  solved  for  CQ0(v)  under  the  previously  stated 
condition  that  only  lower-ordered  feedback  modes  are  significant.  Then 
(21)  can  be  solved  for  I00(v)  • 

2.2  PLANE-WAVE  EXPANSION  MOMENT  METHOD 

j)J) 

The  array  £  and  active  impedances,  z  (se>s^),  z  (s^.s^),  array  £  to 

fib 

array  £  mutual  impedance  z  (s  ,s  ),  and  the  array  £  reference  element  short- 

cx  ^ 

circuit  port  current  Igg,  are  considered  in  this  section.  The  development  is 
not  limited  to  array  radiating  elements  oriented  in  the  plane  of  the  array. 
Scattering  from  element  feedlines  and  radiation  and  scattering  from  elements 
having  arms  inclined  with  respect  to  the  array  plane,  thus,  are  accounted  for. 

Consider  first  array  b_  with  the  ground  screen  absent.  A  port  voltage 
or  current  source  array  excitation  that  is  uniform,  except  for  the  pro¬ 
gressive  phase  -jk  (md  s  +  nd  s  )  imparted  at  the  mn  element  with 

e  e  y  y 

respect  to  the  reference  element  (m  -  n  -  0),  results  in  the  current 


tmn(?,+?«»>  "  V?,)  exp  l'Jk,BdeSe  '^Vy1  (23) 

along  the  nmth  element  wire  axis,  where  r'  is  a  point  along  the  reference  ele¬ 
ment  wire  axis  (path  C'  of  Figure  9)  and  r  is  the  vector  from  the  reference 

element  feed  to  the  mn  element  feed.  This  form  for  t  (r'+r  )  also  will 

mn  mn 

arise  if  the  array  is  excited  by  a  plane  wave  with  directional  cosines 

sy,  3^,  s ^  for  the  propagation  direction.  Under  this  excitation,  once  Iqq(t') 

is  determined,  the  current  on  any  array  element  is  given  by  (23). 
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-  PATH  C 

(WIRE  SURFACE) 


Figure  9.  Curved  Thin  Wire  Array  Reference  Element 
Modeled  b.y  Linear  Segments 


22 


The  boundary  condition  from  which  I^Cr')  is  determined  is 


£*X  -  E 

tan  tan 


along  the  surface  of  the  reference  element,  where  E  is  the  field  radiated  by 
all  I  (r '+r  ) »  E6*  is  an  exciting  plane  wave  in  the  absence  of  the  array 

(impressed  field),  E  is  the  total  field,  and  the  subscript  "tan"  denotes 
"component  tangential  to  wire  surface."  Along  perfectly  conducting  wires, 

-V  -* 

*  0  at  all  points  except  at  a  source  port,  where  E  is  the  source 
tan  r  r  v  '  tan  h 

field  of  an  ideal  port  voltage  source,  V^,  and  it  is  assumed  that  each  mn 
element  is  excited  by  port  voltage 


V  =  V_n  exp  [-jkmd  s  -jknd  s  ] 
mn  00  r  J  eeJ  yy 


A  partial  difference  operator,  point  matching  moment  method  solution 
of  the  thin  wire  approximation  of  (24)  results  in  an  approximation  of 
Ioo^r')  [2].  In  this  procedure,  the  wire  axis  of  the  reference  element  is 
approximated  by  an  Ng  straight  line  segment  fit  (path  C’  of  Figure  9)  begin¬ 
ning  half  a  segment  in  from  wire  ends  where  the  current  is  zero.  The  refer¬ 
ence  element  may  be  composed  of  disjoint  wires  (as  with  parasitics)  and  also 

may  have  multiple  wire  junctions  (as  with  top-loaded  dipoles).  First, 

-►  •> 

IqO^')  is  expanded  as 


.-r,  . 

Ioo(c  } 


h  Vr'>h 


where  fj(r' )  is  a  unit  "pulse"  of  wire  current  traversing  the  jth  segment  of 

length  A ^  ,  2.^  is  the  j6*1  segment  of  path  C'  reference  direction  unit  vector 
and  Ij  is  an  unknown  coefficient  to  be  determined.  The  substitution  of  (25) 
into  (24)  results  in 


Eh 
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'*55  -►  t  h 

vihiite  E  (t,j)  is  the  field  radiated  by  the  infinite  array  ot  j  segment  expan¬ 
sion  functions 

f.(r')£.  exp  f-jkmd  s  -iknd  s  ] 

33  eeJ  yyJ 

A  matrix  equation,  solvable  for  the  I  ,  results  from  (26)  by  a  point  matching 
procedure.  Equation  (26)  is  satisfied  at  the  centers  of  each  segment  along  a 
path  C  parallel  to  the  segmented  path  C'.  This  is  equivalently  an  impulse 
weighting  moment  method  whereby  the  weighted  average  of  (26)  along  C  is  taken 
with  respect  to  each  of  the  weights 

w .  (r )  =  A.  6(r-r.)i. 

i  l  ii 


where  5(r-r.)  is  the  unit  impulse  function  and  r.  is  the  position  vector  of 

t  h  * 

the  center  of  the  i  segment  of  path  C.  The  result  is 


(7.11  =  Vex  +  V 


(27) 


th  —  ”*ex  _ 

where  the  i  elements  of  the  N  x  1  column  vectors  T,  V  ,  and  V  are, 

s  ’  ’ 

respectively,  1^ 


1. 

i 


/-+ex  -+  -*■  ->ex  -> 

E  -w.  (r)  di.  =  A .  E  (r.)  ■%. 
l  ill 

c 

->  ->  ~ 

E-w.  (r)  dt  =  A.  E(r . ) -l . 

l  iii 


(28) 


and  the  i j ( element  of  the  M  x  N  nome-i.t  matrix  [7.1  is 

s  s 


il 


/  *  ->  s  ► 

w.(r)-E  ( r , j 


)  di  =  -  A  .  ES  { r  .  ,  j )  •  i  . 

i  l  i 


(29) 
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In  determining  Z  the  field  E^Cr^.j)  is  represented-  in  terms  of  a  vector 
potential  A  and  gradient  of  a  scalar  potential  <t> .  Following  the  development 
in  [2,  Chapter  4],  the  scalar  potential  term  is  approximated  by  a  partial  dif¬ 
ference  for  the  gradient  resulting  in 


Z  =  juA  l  *A(r  , j)  +  <t(r  ,j)  -  $(r  ,j) 


(30) 


where 


r.  =  r.  +  (A./2U. 


iii 


r.  =  r.  -  (  A./2U. 
1,  1  1  1 
D 


and  to  is  the  radian  frequency.  The  vector  potential  arises  from  the  infinite 
array  of  progressively  phased  unit  amplitude  pulse  segments 


fj(r')£j  exp  [-jkmdese  -jknd^s^] 


and  the  scalar  potentials  from  partial  difference  approximations  to  the  linear 

charge  densities  associated  with  the  infinite  array  of  f  (r’ ).  The  reference 

->■  -►  -J 

element  charge  density  Oj(r')  associated  with  fj(r')  becomes  two  "shifted" 
pulses  of  charge  density,  each  of  length  A  ^ .  Thus 


o-(r’) 


f.(r') 

■^a 


f.  (r') 


juA . 
3 


(31) 


where  fj  is  f^  shifted  in  the  direction  to  center  at  rj  -  r j  +  (Aj/2)ij 

S  3l 

(rj  is  the  position  vector  of  the  center  of  the  j  segment  of  path  C')  and 

f  is  f,  shifted  in  the  direction  -i  to  center  at  r'  =  r!  -  (A  /2)i  . 

J  |j  J  J  J  J  J  J 

Plane  wave  expansions  for  £(r  .1)  and  $(r  ,i)  are  derived  below. 

a  a 

Munk  and  Burrell  obtained  the  expression  [1,  Eq.  3]  for  the  field  from 
an  infinite  planar  array  of  point  current  dipoles.  They  indicated  that  this 
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expression  is  valid  tor  arbitrarily  oriented  dipoles;  tiie  dipoles  need  not  be 
oriented  parallel  to  the  array  plane.  Munk  and  Burrell  also  mentioned 
that  a  moment  method  can  be  developed  from  this  expression.  However,  the 
mutual  impedance  expression  they  derived  (1(  ._q,  p 9 ]  restricts  mutually  cou¬ 
pled  antenna  segments  to  have  disjoint  z  axis  projections.  The  generalization 
of  the  moment  method  to  handle  "out  of  plane"  elements  necessitates  more  com¬ 
plicated  expressions.  these  are  derived  here. 

Consider  two-dimensional  arrays  of  incremental  current  dipoles  and  incre¬ 
mental  linear  charge  densities  having  ,wt*'  elements  given  respectively  by 

f  (r'+r  )d*’  -  f .  (r’).U'  v  .expf-jkmd  -  -  j had  s  j 

mn  mn  j  j  ■  «  c  v  y 


o  (r'+r 
tnn  mn 


)  d  ’ ' 


,  (r  '  ;  J '  exp  [  -  j  kmd  s  --  jtnd  *-  : 


Following  the  work  of  Munk  and  Burrell  [1,  Appendix  A  1  and  Kornbau  [4],  the 
corresponding  vector  potential  dA(v  ,j)  and  scalar  potential  d-(r^,j)  at  a 
point  r  are  given  by  the  plane  wave  expansions 


dA(ra,j) 


dl'n  f  .  <r ' )  ;  . 

_ J  J 

2ju)d  d 

^  y 


E  E 


rj=_co 


exp  [ -jk(r  -r ’ ) * g^j 


(32) 


df’o  (r') 

d'(yJ)  =  2iki -r 

x  * 


E  E  - 


xp  f - j k ( r  -r’ ) -g+] 


(33) 


for  (r  -r')-z  ? 0  where 

n 

n  =  free  space  wave  impedance 
=  free  space  permitivity 

8-  "  xfVr'  d  "  d'd  •  1  +  yfsv+fi  i  ]iz*z 

x  xv  y 


y~i~  t  s 


\  '  ;.v  2  t  •*.  ? 

+p  j- •  +n  - — T'  ]  ~  I  ^  +Q  /  *1 

x  cl  1  J  d  v  o 

x  xv  V 


i  =  free  space  wavelength 


~k.  I — 
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If  the  expression  under  the  radical  in  gz  is  negative,  g^  is  chosen  to 
be  negative  imaginary;  if  the  expression  is  positive,  gz  is  nonnegative  real. 
This  assures  that  all  plane  wave  modal  fields  travel  outward  from  the  source 
array  plane  and  nonpropagating  modes  attenuate.  A  plus  sign  subscript  on  g 
corresponds  to  field  points  r^  in  front  of  the  array  plane;  i.e.,  (r^r' ) *z  > 0. 
A  minus  sign  subscript  corresponds  to  field  points  behind  the  array  plane; 
i.e.,  <0. 

Since  fj(r')  is  the  reference  element  segment  unit  "pulse"  and 

r*  =  r (as  defined  in  Figure  10),  (31),  (32),  and  (33)  yield,  after 
interchange  of  integrations  and  summations, 


A(ra,j)  = 


f 

C' 


dA (r  ,  j ) 


i . nA . 

— 3-JL_  V 

2jtod  d  /  a 
xy 

p=-“ 
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The  iJj  (a,  6)  functions  assume  one  of  three  different  forms: 

pq 


(33) 


exp[-jk(r  -r' )  -g  )]  sin[k  (X. .  *g  )  A  ,/2] 

i i  (a,  6)  =  - 2 — 2 - +  _  -  i  +2 - 

P3!  g  k(d  -g  )A  /2 


(36a) 


it  (vip-iifISj-i! 


exp[-jk(r  -r' ) -g  ]  sinfk(£.-g  )A./2] 

♦  ■  (a,  6)  =  - 2—2 - - - 3.—  -J - 
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Figure  10.  Expansion  and  Weighting  Seaments 


Equation  (36a)  applies  when  the  entire  source  segment  is  closer  to  the 
ground  screen  than  is  the  field  point  The  plus  sign  subscript  is  indica¬ 

ted  since  the  plane  waves  then  are  propagating  in  the  positive  z  direction. 

# 

Equation  (36b)  applies  when  the  entire  source  segment  is  further  from  the 
ground  screen  than  is  r^.  The  negative  sign  subscript  now  appears  since  the 
plane  waves  are  propagating  in  the  negative  z  direction.  Equation  (36c) 
applies  when  r^  defines  a  point  common  to  a  plane  that  simultaneously  di¬ 
vides  the  source  segment  and  is  parallel  to  the  array  plane.  The  i  coor¬ 
dinate  along  the  source  segment  of  this  point  is  defined  by  equating  z  axis 

projections  of  r  and  r'  +  Vi ,  yielding  l  given  by  (37)  (Figure  10  with 

•*  ■+  ®  “  3 

ra  »  r^  )  .  For  ^**<0,  e.g.,  the  direction  indicated  in  Figure  10, 
integration  with  respect  to from  -A^ /2  to  1  involves  plane  waves  propaga¬ 
ting  in  the  negative  z  direction;  hence,  the  -  subscript  is  indicated  for  9 
in  the  second  set  of  terms  in  (36c).  For  the  integration  from  1  to  Aj/2,  the 
properly  chosen  plane  waves  propagate  in  the  positive  z  direction;  hence,  the 
first  set  of  terms  in  (36c).  Equation  (36c)  also  applies  when  £j*z>0.  The 
upper  (lower)  of  the  ±  signs  in  (36c)  applies  for  ij’z>0  (£^*z<0).  If 
£  -z  -  0  (I  not  defined),  either  (36a)  or  (36b)  applies. 
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Equations  (30),  (34),  and  (35)  yield 


A . A . n 

,  =  _LJ _ 

13  2d  d 

x  y 
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EE 
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v  ci.j) 
pq 
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EE 
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pq  a  a  pq  a  b 


1I1  (i:  ,  j  )  +  ii<  (i,  ,  j.  ) 
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In  (38),  the  field  point  vectors  r.,  r.  ,  and  r.  occur  (not>;-  defining  equa- 

1  La  Lb 

tions  following  (30)).  These  vectors  define  the.  midpoint  and  endpoints  of 
the  itn  weighting  segment  along  path  C  as  shown  in  Figure  10.  Thus 


>  -►  -*■ 

r  .  =  r '  +  a  . 
111 


where  ja.'|  is  the  wire  radius  a.  along  the  i  segment.  The  vector  direction 

-*■  .  4 

of  a.  is,  of  course,  normal  to  l: .  The  vector  a.  is  defined  fully  by  requir- 
1  *  1  -  1 
ing  that  a^  also  lie  in  the  plane  formed  by  and  z.  This  speeds  convergence 

of  the  infinite  summations  in  (38)  for  segments  oriented  nearly  parallel  to 

the  xy  plane.  This  condition  is  satisfied  for 


at  =  a^  x  (z  x  x  '(2  x  O  ! 


If  a  ground  screen  is  present  a  distance  d  behind  the  array,  image 
fields  must  be  considered  In  each  impedance  element  computation.  Let  double 
primed  position  vectors  and  segment  path  lengths  refer  to  those  for  image  cur¬ 
rents  and  associated  charge  densities.  Then,  since  the  coordinate  origin  lies 
on  the  ground  screen. 


r =  x(x-r!)  +  y(yr!)  -  z(z-r!) 
J  J  1  J 


f-V  =  x(x-!.j)  +  y(y •  t  )  -  z(z*iL) 


t  h 

and  the  ij  element  of  the  moment  matrix  for  an  infinite  array  over  a  per¬ 
fect  ground  screen  becomes 


ZgS  -  z  Z,lmage 
ij  iJ  ij 


(39) 


where  Z  is  given  by  (38)  and  Z^^^is  7  of  (38)  with  (a,B)  of  (38)  com- 

ij  ij  iJ  pq 

^  ^  A  A 

puted  by  (36)  with  r^  replaced  with  rV,  and  iL  replaced  with  2. V . 

For  plane  wave  incidence,  the  voltage  column  vector  elements  given 
by  (28)  must  be  appropriately  modified  to  account  for  reflection.  If 

EfiX  -  EeX  i  +  EeX  y  +  EeX  £ 
x  y  J  z 

6X  6X  SX 

(where  E^  ,  E^  ,  and  E^  are,  by  the  plane  wave  condition,  constrained  by 

GX  GX  6X  til 

sE  +sE  +sE  =0),  the  i  element  of  the  voltage  column  vector 
xxyyzz 

including  reflection  from  the  perfect  ground  screen  is  given  by 


-ik(s  x+s  y) 

v;x'g^Vi-e  x  y 


-2jsin(ks  z)(Eex  x+EeX  y)  +  2cos(ks  z)Eex  z 
J  z  x  y  J  z  z 


(40) 


Equation  (27)  then  becomes 


[ZgS]I  -  vex’gs  +  V 


(41) 


Equations  (36),  (38)  ,  and  (39)  simplify  if  the  radiating  elements  lie  in 
the  array  plane.  For  this  special  case,  (36)  and  (38)  reduce  to 
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"ij  2d 


A  2-rZ-r(AiVVV 


p=-oo 


(42) 
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Only  the  "+”  sign  for  g+  appears  in  (42)  and  (43)  because  of  the  simplifying 
assumption  of  wires  lying  parallel  to  the  xy  plane  and  the  path  C  is  chosen 
to  lie  a  wire  radius  in  front  (positive  z  direction)  of  C'. 

The  progressive  phasings,  s  and  s  ,  which,  for  a  particular  p  and  q, 

x  y 

cause  g  to  be  zero,  can  result  in  an  infinite  Z...  These  progressive  phas- 
z  ij 

ings  correspond  to  "grating  lobe  singularities"  since  they  correspond  to 
grating  lobes  just  entering  the  visible  range;  i.e.,  nonattenuating  plane 
waves  propagating  along  the  array  plane.  This  is  one  possible  source  of 
"blind  spots"  and  rapidly  varying  active  impedances  with  scan.  This  is  only 
a  "potential"  condition  for  infinite  Z„  since,  although  is  in  the  denom¬ 
inator  of  '!>  ,  in  some  situations  a  compensating  effect  occurs.  Such  is  th 

pq  . 

case,  for  example,  in  E-plane  scans  with  periodic  infinite  arrays  of  v  di¬ 
rected  dipoles  where  E-plane  grating  lobe  singularities  are  suppressed.  This 
is  clear  from  (42)  for  small  k.A^/2  since,  in  the  E-plane,  g+-x  =  0,  I^-g+  =  1 
as  g^->-0,  and  the  parenthetical  term  multiplying  ij^(i,j)  vanishes. 

The  presence  of  a  perfect  ground  screen  also  inhibits  blind  spots  due 
to  grating  lobe  singularities  if  the  radiating  elements  lie  parallel  to  the 
array  plane.  If  d  is  the  height  of  the  array  above  the  ground  screen  and 
the  positive  z  axis  is  directed  from  the  ground  screen  to  the  array  as  for 
array  b,  7.^?  of  (39)  is  given  by  Z„  of  (42)  with  ii^(i,j)  given  by 


pq 


(i,j)  = 


2j  sin (kdg  ) 


exp  [- jk  (  (r. -r.' ) -g  +dg  )  ] 
i  j  +  z 


sin[k (£. -g  )A  /2] 

- --4 — — J -  (44) 

Mi.-g+)4./2 


in  place  of  (43).  It  is  apparent  from  (44)  that  is  bounded  at  a  grating 

pq 


sin  kdg 

lobe  singularity  since  -  is  bounded  as  g  4  0.  This  ground  screen 

^z  ^ 

suppression  of  a  blind  spot  at  a  grating  lobe  singularity  assumes  the  radia¬ 
ting  elements  lie  in  the  array  plane.  If,  for  example,  the  dipole  arms  are 
inc Lined  with  respect  to  the  array  plane  or  feedline  scattering  is  present, 
the  ground  screen  may  enhance  the  blind  spot  effect  at  a  grating  lobe  singu¬ 
larity  above  what  it  would  be  if  the  ground  screen  is  absent;  thus,  consider¬ 
ation  of  out-of-plane  radiating  elements  via  the  general  expressions  (36), 
(38),  and  (39)  is  important. 
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The  array  active  impedance,  z'3^3(s 

ing  (41)  for  the  source  segment,  i  =  i  , 

VjX,gs  ■  0  for  all  i  and  V.  =  0  for  all 
i  i 


e,s  ),  then  can  be  determined  by  solv- 

current  coefficient,  1^  ,  with 

i  4  i  .  Then  s 

s 


bb 

z 


(s„ 


The  array  a^  active  impedance,  zaa(se>s^),  can  be  obtained  in  the  same  manner. 

The  array  a  reference  element  short-circuit  current,  Iqq.  can  be  deter¬ 
mined  by  solving  (41)  for  I 


ments  of  vex’®s  are  given  byS(40). 


Now,  however,  V  =  0  for  all  i  and  the  ele- 


ab. 


The  array  b^  to  array  £  ground  screen  penetration  coupling,  z  (s^.s  ), 
feet  ground  screen  can  be  determined  from  (41)  where  [Z®S]  is  the  array  a 
ground  screen  present  moment  matrix,  V  =  V  +  V,and  the  elements  of  V 

“►t 

are  given  by  (28)  with  E  replaced  by  the  field  E  transmitted  through  the 
ground  screen  with  array  b^  excited  by  current  sources  given  by  (4).  (A  means 
for  determining  E C  is  discussed  below).  The  element  in  the  solution 

vector,  I  of  (41),  is  then  the  short-circuit  array  a  reference  element  port 
current.  The  reference  element  open-circuit  voltage  is 


aa ,  \ 

-z  <se,sy) 


thus. 


ab 
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V 


aa/  ■,  ra 

z  (se>V  Xi 


00 


(45) 


The  field  is  obtained  readily  from  knowledge  of  the  plane  wave 
transmission  coefficients  of  the  ground  screen  and  the  plane  wave  decom- 
postition  of  the  field  from  array  1).  The  array  b_  current  is  determined  with 
array  l>  excited  with  port  current  sources  (4).  The  array  b  reference  element 
current  vector  I*5  is  the  solution  I  to  (41)  with  V®x,gs  =  0  for  all  i  and 
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h_  =  (it.  X  g.)  x  g_ 

can  be  derived  from  [1,  Eq.  3  ).  The  subscripted  negative  sign-,  indicate  that 

(46)  applies  for  (r  -r ! )  •  z  ■  0  as  is  appropriate  since  the  ground  screen  is 

->i  -►t 

"behind"  array  b  .  Since  F.  is  a  superposition  of  plane  waves,  E  is  obtained 

by  applying  the  proper  transmission  coefficients  to  each  component  plane  wave. 

Each  plane  wave  in  (46)  propagates  in  the  g  direction;  thus,  g_  and  z  (ground 

screen  normal)  form  the  corresponding  "plane  of  incidence."  A  unit  vector 

normal  to  the  incidence  plane  is 

g  x  z 
!g_  x  z| 

an-!  a  unit  vector  parallel  to  the  incidence  plane  and  normal  to  g  is 

n  ||  =  g_  x  nj_ 
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Any  plane  wave  travelling  in  the  g_  direction  can  be  decomposed  into  two 

-N  -> 

plane  waves,  one  with  an  E  field  polarized  along  nj^  and  the  other  with  an  E 
field  polarized  along  n||  •  If  Tj_  and  T  |j  are  the  corresponding  transmission 
coefficients, 


E  (r)  = 


S  OC  90 

j  =  l  J  p=-oo  q=-oo 


niTj^(nj_-h_) 


+  nll  Tll  (nll  -h-1  ^pq 


The  substitution  of  (47)  for  E  in  (28)  yields 
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P=— 00  q=-a 


and  z  (s  ,s  )  is  determined  as  previously  discussed  in  the  paragraph  con- 

^  y  ^  k 

taining  (45).  In  (48)  ^^(i.j)  is  ^  ( j )  with  r  ^  replacing  r. 
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SECTION  I 


IMPLEMENTATION 

The  computer  programming  to  date  of  the  lens  simulator  is  discussed  be¬ 
low.  Two  main  programs  have  been  developed  -  FAM  and  HAM.  The  first  approx¬ 
imation  method  is  carried  out  in  FAM  (Section  3.1).  The  higher  approximation 
method  is  carried  out  in  HAM  (Section  3.2).  The  program  HAM  relies  on  data 
generated  by  FAM  that  is  stored  in  disk  files.  Several  subprograms, 
discussed  in  Section  3.3,  are  used  by  both  FAM  and  HAM.  The  discussions  and 
flow  diagrams  follow  closely  the  theory  nomenclature  of  Section  2  and  [5). 
Computational  considerations  of  several  equations  and  convergence  rational 
are  discussed. 

User  data  input  is  described  in  Section  3.3.11.  Data  is  inputted  in 
an  interactive  manner;  the  computer  queries  the  user. 

The  general  scattering  parameter  model  derived  in  Section  2.1  is  not 
yet  implemented.  This  implementation  is  presently  taking  place.  The  spe¬ 
cialized  scattering  parameter  model  that  applies  to  passive  phase  shifters 
composed  simply  of  lengths  of  transmission  line  is  implemented.  Also  a 
specialized  impedance  parameter  model  [5,  Section  3.2.1]  is  implemented. 

The  user  selects  the  model  most  suitable. 

An  efficient  folded  dipole  radiator  model,  not  discussed  in  the  theoret¬ 
ical  developments,  is  implemented  as  a  user  option.  This  model  is  described 
in  Section  3.1. 

Computer  run  time  is  consumed  primarily  in  the  computation  of  the  infi¬ 
nite  array  moment  matrices.  Atlantic  Research  Corporation  is  investigating 
eliminating  this  computation  from  FAM.  Instead ,  a  third  computer  code  would  be 
constructed  that  generates  the  infinite  array  moment  matrices  for  several 
scan  angles.  The  program  FAM,  then,  would  read  these  matrices  from  data 
files  and  interpolate  to  arrive  at  the  moment  matrix  for  a  specified  scan 
angle.  The  moment  matrix  generating  code  would  best  be  run  in  a  batch-like 
mode  and  FAM  and  HAM  in  a  user- interact ive  mode. 


3.1  First  Approximation  Method 


The  first  approximation  method  (FAM)  implementation  is  discussed  here. 

The  radiating  element  currents  and  port  impedances  are  computed  for  each  cell 
assuming  infinite  periodic  lattices.  The  first  approximation  lens  radiation 
patterns  are  computed  from  these  currents.  The  underlying  theory  is  described 
in  Section  2  and  [5]. 

The  subroutines  called  by  FAM  are  listed  in  order  of  first  invocation 
in  Figure  11.  Subroutines  called  by  subroutines  of  FAM  also  are  listed 
(indented)  beneath  their  respective  calling  programs.  A  flow  diagram  of 
FAM  is  given  in  Figure  12 . 

Lens  data  are  read  in  via  LENSIN,  cell  data  via  CELLIN  and  cell  illumi¬ 
nation  via  FEED.  Subroutine  FEED  assumes  a  point  source  located  a  global  z 
axis  distance  behind  the  lens.  The  effective  radiated  power  and  E-field  x,y 
polarization  components  (global  system)  at  the  cell  coordinate  origin  must  be 
provided  for  each  plane  wave  component  of  illumination.  (The  illumination  is 
the  field  in  the  absence  of  the  lens;  the  global  coordinate  origin  is  in  the 
lens  ground  screen) . 

An  option  for  modeling  folded  dipole  radiating  elements  efficiently  is 
available  in  FAM.  If  the  option  is  exercised,  the  user  supplies  the  folded 
dipole  transmission  line  mode  characteristic  impedance  ZQ  and  dipole  length 
(end  to  end)  L.  The  transmission  line  mode  impedance  is  computed  in  FAM  as 

Zt  -  jZQ  tan(~  kL) 


This  impedance  is  combined  with  the  infinite  array  folded  dipole  radiation 
mode  impedance  z  to  yield  the  folded  dipole  scan  impedance  given  by  [6] 


folded  dipole  scan  impedance 


4Ztz 

Zt+2z 


The  infinite  array  folded  dipole  radiation  mode  impedance  z  is  the  scan  imped 
ance,  including  scattering  from  feedlines,  of  an  array  of  simple  dipoles  of 
equivalent  folded  dipole  radii  The  wire  geometry  supplied  by  the  user 

that  specifies  the  radiating  element  must,  therefore,  treat  the  folded  dipole 
as  a  simple  dipole  of  radius  ag.  A  good  approximation  for  ag  can  be  obtained 
through  two  dimensional  analyses  as  indicated  by  Wolff  [6,  pg.  61].  A  first 
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Figure  11.  Program  FAM  Subprograms 
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CALCULATE  NORMALIZATION  CONSTANTS 
FOR  MODULE  S-PARAMETERS. 


CALCULATE  PHASE  DELAY  ANO  CHARACTERISTIC  IMPEDANCE 
OF  REFERENCE  ELEMENT  PASSIVE  LINE  PHASE  SHIFTER 


BEGIN  LOOP  OVER  FEEDBACK 
MODES  Ua.  Uh 


Figure  12,  Flow  Di 
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for  FAM 


N» 


TRANSMISSION  LINE  MODULE  PHASE  SHIFTERS 
ASSUMED.  CAL CUt ATE  ARRAY  a  AND  ARRAY  b 
REFERENCE  ELEMENT  ZERO  FEEDBACK  MODE 
PORT  CURRENTS  lg0(0),  1§0(0)  AND 
SCATTERING  VARIABLES  Cg()(()), 


IS  GROUND  SCREEN  FEEDBACK  PRESENT? 


YES 


MUTUAL 


CALCULATE  GROUND  SCREEN  PENETRATION  MUTUAL 
IMPEDANCE  COUP.  TNG  zab  BETWEEN  ARRAY  b 
AND  ARRAY  a.  AlSO  COMPUTE  ARRAY  a  EXCI¬ 
TATION  VOLTAGE  AND  SHORT  CIRCUIT  CURRENT 
COLUMN  VECTORS  DUE  TO  CURRENTS  ON  ARRAY  b 


IF  MODULES  (PLUS  FEEDLINES)  REPRESENTED 
BY  TWO  PORT  IMPEDANCE  PARAMETERS,  FORM 


iab 


•  C  -  >“ 


CALCUI  ATE  FEEDBACK  MODE 
CONTRIBUTION  TO  ARRAY  a 
AND  ARRAY  b  REFERENCE 
I U Mi  NT  PORT  CURRENTS 


(  FIRST  FEEDBACK  LOOP?  - (5) 


NO  I 

Q 


Figure  12.  Continued 
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Figure  12.  Concluded 
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approximation  to  a^  for  a  folded  strip  dipole  with  parallel  strips  separated 
a  small  distance  with  respect  to  the  overall  folded  dipole  width  w  is 

a  =  w/4 
e 

The  user  may  choose  either  an  impedance  parameter  model  or  a  scattering 
parameter  model  for  the  module  (plus  feedlines)  two-port  implementation.  The 
impedance  parameter  model  is  based  on  the  derivation  in  [5,  Section  3.2.1]. 

It  allows  for  active  modules  (amplifiers)  but  restricts  the  impedance  param¬ 
eters  to  an  idealized  model;  impedance  parameters  must  be  uniform  functions 
of  progressive  phase  angle  setting  except  for  the  phase  of  the  forward  mutual 
impedance  parameter.  That  phase  varies  from  modple  to  module  by  a  progres¬ 
sive  phase  increment.  Be  advised  that  the  impedance  parameter  model  de¬ 
scribed  in  [5]  is  based  on  [5,  Eq.  (3-12)]  whereas  the  development  in  this 
report  is  based  on  (6).  The  two  equations  differ  by  a  feedback  mode  offset, 
and  the  restriction  to  only  nonnegative  feedback  modes  in  [5].  The  offset  im¬ 
plies  that  the  u  *  p  mode  computed  using  (6)  is  the  u  *  p-1  mode  computed 
using  [5,  Eq.  (3-12)].  The  mode  incrementation  loop  in  FAM  follows  the  (6) 
convention.  Thus,  if  Impedance  parameters  are  chosen,  the  u  *  p  mode  compu¬ 
tation  of  Ib  in  FAM  should  be  interpreted  according  to  [5,  Eq.  (3-12)]  as 
mn 

the  u  =  p+1  mode. 

If  Mutual  is  bypassed  (ground  screen  assumed  perfectly  conducting)  only 
a  scattering  parameter  module  representation  is  allowed. 

The  scattering  parameter  representation  for  the  module  plus  feedlines 
is  presently  implemented  for  the  special  case  of  modules  composed  simply  of 
lengths  of  transmission  line.  The  characteristic  impedances  of  the  module 
transmission  line  and  feedlines  are  assumed  equal  (*‘Rg).  Equations  (21)  and 
(22)  then  apply  for  determining  the  array  a_  and  array  b^  reference  element 
port  scattering  variables  Cqq^  and  currents  for  "feedback"  modes 

u  =  0,  ±1,  ±2 ,  . . .  .  Equations  (22)  generally  are  not  decoupled.  A  simpli¬ 
fication,  employed  in  the  present  implementation,  is  that  only  nonnegative 
feedback  modes  u  are  significant.  Equations  (22)  then  decouple  to  yield  the 
recursion  relations 
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=00(0)  •  (/ro[°l 


500(u)  •  «ol">  +  '*J>  «*u>  -  V”»  IW!500(ti-l> 


U  -  1,  2, 


The  corresponding  port  currents  are  found  from 


j  _ L_  c 

00(0)  r~  L00(0) 


*00(u)  =  "  y—  (G00(u)  "  [W]C00(u-l))  u  1*  2’ 


A  qualitative  justification  for  omitting  negative  feedback  modes  is  given  in 
[5(pg.  3-4)].  The  validity  of  this  assumption  lies  in  observing  convergence 
in  Iqq(u)  with  increasing  u. 

The  above  four  equations  are  implemented  by  first  computing 
_  _  .  aa  _  .  .  ba 

Q11  =  Gll(zu  “  V  G12  Zu 

Q  *=  G  zab  +  G  ,(zbb  -  R  ) 

12  11  u  12  u  0 

#  aa  v  ba 

Q21  "  G21(zu  "  R0)  +  G22  Zu 


Q22  *  G21  Zu  +  G22(zu  "  V 


,  aa  aar  ,  _  .  _  bb  ab  ,  ba 

where  z  ■  z  [s  +ua  ,  s  -Hia  ]  and  so  forth  for  z  ,  z  ,  and  z  ; 
u  e  e  y  y  u 


Ro  +  zbb 

,  _  _ 0  u _ 

'll  .  bb.  aa.  ab  ba 

(R.+z  )  (R_+z  )  -  z  z 

0  u  0  u  u  u 


12  ,  aa.  bb.  ab  ba 

(R_+z  )  (Rn+z  )  -  z  z 

0  u  0  u  u  u 
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ba 


21  /n  ,  aa,  J  bbN  ab  ba 
(R.+2  )  (R-+z  )  -  z  z 

0  u  0  u  u  u 


G22 


R  +  zaa 
0  u 


(Rn+zaa)  (Rn+zbb)  -  z3b  zba 
0  u  0  u  u  u 


and  provision  for  a  feed  "forward"  imperfect  ground  screen  penetration  trans¬ 
fer  impedance  zba,  in  addition  to  the  feedback  term  zba,  is  included.  It  is 
ba 

expected  that  zy  will  he  set  to  zero  during  an  actual  lens  analysis  for  the 
reasons  cited  in  Section  2.1  concerning  (1)  and  (2).  During  validation  of 
the  simulation  program,  however,  provision  for  :ba  other  than  zero  m?.\  oe 
helpful.  Then 
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and,  for  u  =  1,  2,...  , 
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where  i.  is  the  length  of  reference  element  transmission  line  (module  plus 
feedlines). 

If  the  ground  screen  is  perfectly  conducting, 

\  "  Zu3  "  G12  "  G21  "  Ql2  =  *21  "  °  and 

I00(u)  =  ° 

T00(u)  =  °  U  °dd 
3.2  Higher  Approximations  Method 

The  theory  underlying  the  higher  approximations  method  (HAM)  is  given 
in  [5,  Section  3.3].  First  approximation  radiating  el emeu l  currents 
initialize  an  iterative  procedure  that  updates  the  current  distributions 
on  elements  near  discontinuities  in  periodicity  (lens  edge,  sec.  ion  interfact 
failed  modules,  etc.)  The  lens  radiation  pattern  also  is  updated. 

The  present  implementation  of  HAM,  discussed  here,  computes  only  the 
second  approximation.  This  is  expected  to  yield  accurate  results  for  lattice 
spacings  exceeding  a,  0.251. 

The  subroutines  called  by  HAM  are  listed  in  order  of  first  invocation 
in  Figure  13.  Subroutines  called  by  subroutines  of  HAM  also  are  listed 
(  indented  )  beneath  their  respective  calling  programs.  A  flow  diagram  ol 
HAM  is  given  in  Figure  14. 

The  combined  module  and  feedlines  (transmission  line  mode)  two-port 

model  can  be  specified  in  terms  of  the  specialized  impedance  pat  3  used 

in  FAM  or  general  scattering  parameters.  When  impedance  parameters  art 

(2)  a 

specified,  the  array  a^  and  array  second  approximation  port  currents  1^" 

Iq2 ) b  for  a  0^  array  element  are  computed  by  solving  [5,  '~qs.  '3  ■p/p  n 

When  scattering  parameters  are  specified. 
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Figure  13.  Program  HAM  Subprograms 
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GEOMET 


Figure  14.  Flow  Diagram  for  HAM 
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LOOP  OVER  COLUMNS,  na,  OF  ARRAY 

(na  INCREASES  ALONG  THE  e  LATTICE  COORDINATE) 


CALCULATE  iLn  ELEMENT  xy  COORDINATES  WITH  RESPECT  TO  CELL  COORDINATE  ORIGIN 


xi  =  mdx  +  xaref  *i  =  ndy  '  +  ~va 


CALCULATE  i tn  ELEMENT  LOCATION  IN  GLOBAL  COORDINATES 

x.  =  x  +  x.  cos  t  -  y.  sin  $ 

ig  cref  i  'q  m  'g 

y.  =  y  ,  +  x.  sin  <■  +  y.  cos  a 

-Mg  Jcref  i  g  vg 


LOOP  OVER  CELLS 


INPUT  ELEMENT  PLACEMENT  DATA  FOR  r"  CELL 


LOOP  OVER  ARRAYS  WITHIN  CELL 


DEFINE  REGION  OF  INFLUENCE  OF  RADIATING  ELEMENTS  IN  l 
ARRAY  DESCRIPTION  RELATIVE  TO  i th  RADIATING  ELEMENT 


ARE  ANY  ELEMENTS  WITHIN  REGION  OF  INFLUENCE? 
•  m0  =  n0  =  0? 


INPUT  RANDOM  ACCESS  FILE  RECORD  IDENTIFIERS 
FOR  FEEDBACK  MODE  INFORMATION  OF  CELL 


LOOP  OVER  "COLUMNS"  IN  REGION 
OF  INFLUENCE,  m  =  m.(  ,m2 


LOOP  OVER  "ROWS"  IN  REGION 
OF  INFLUENCE,  n  =  n-,  ,n2 


Figure  14.  Continued 
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Figure  14.  Continued 
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Figure  14.  Continued 
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UPDATE  FIELD  PATTERN 


r 


is  given  by 


I0  -  -  [R]^  «D]  -  ts03)  c0 


where  [R]  and  [U]  are  defined  in  the  discussion  following  (11),  [Sq]  is  the 

zero1"*1  element's  scattering  parameter  matrix,  and  is  the  zero1"*1  element's  in¬ 
cident  scattering  variable  column  vector.  The  vector  C_  is  the  solution  to 


[R]2  <[UJ  +  [Sj  +  [z]  (|U]  -  [Sn]))  Cn  =  [z]  If 


za  0 


0  z 1 


JoC  = 


ex  (l)aa'  T(l)ab 
0  0 


(l)bb’ 


za  and  z  are  the  array  <1  and  the  array  b_  isolated-element-above-a-ground- 
screen  input  impedances;  and  1^,1^  ,  Iq  ,  and  Iq  are  the  zero  ele¬ 

ment's  short  circuit  port  currents  due  to  the  illumination  and  first  approx¬ 
imation  currents  on  neighboring  elements  [5,  Section  3.3.1]. 

3. 3  Subprograms 

The  subroutine  and  function  subprograms,  typically  not  provided  by  a 
computer  system,  that  are  used  in  FAM  and  HAM  are  described  below.  The 
discussion  follows  closely  the  theoretical  development  of  Section  2 
and  [ 5 ] . 
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3.3.1  ZMATG 


Subroutine  ZMATG  fills  the  moment  matrices  [Z]  and  [Z®s]  for  planar 
arrays  of  general  "out  of  plane"  thin  wire  radiators  with  and  without  a  per¬ 
fect  ground  screen.  The  ground  screen  is  located  in  the  xy  plane  parallel  to 
the  array  plane.  Block  Toeplitz  symmetry  allows  reducing  the  number  of  com¬ 
puted  matrix  elements  corresponding  to  interaction  between  segments  on  the 
same  straight  wire  to  the  number  of  segments  on  that  wire.  (Toeplitz  symme¬ 
try  means  Z ^  =  Z^+^  for  p  an  integer.)  Toeplitz  symmetry  generally  is 

not  exhibited  by  matrix  elements  corresponding  to  interaction  between  seg¬ 
ments  on  different  wires.  Also  Toeplitz  symmetry  generally  is  not  present 
in  the  ground  screen  image  interaction  matrix  [Z*1113®]  that  subtracts  from 
[Z]  to  yield  [Z®S].  However,  significant  computational  savings  usually  is 
achieved  since  closely  spaced  segments  usually  require  the  most  computation 
due  to  slowly  converging  plane  wave  expansions  and  closely  spaced  segments 
usually  are  located  on  the  same  wire. 

The  ij1"*'  element  of  the  free  space  (ground  screen  absent)  moment  ma¬ 
trix  [Z]  is  computed  in  OFFD  for  i  4  j  (off-diagonal  element)  as  defined  by 
(38).  There,  PSIGEN  is  invoked  five  times  to  approximate  the  five  plane  wave 
expansions 


iKa,  6) 


00  CO 


EE 


’J'  (a»  6) 

pq 


(49) 


qs— oo  p=— oo 


where  a, 6  =  i,j;  ifl , j a ;  ±a, jfe ;  ife , j a ;  or  ib,jfe.  The  i,jth  element  of  [ Z) 
is  computed  there  according  to 


Z.  . 


ij  2d  d 
x  y 


Vj(VV  *(i’J)  --T  *(ib’V 

J  J  V  V 


+  7  *(ia’jb)+7  *(W 


(50) 
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Equation  (38)  simplifies  for  diagonal  elements  of  [z].  The  i  diag¬ 
onal  element  is  computed  in  ZMATG  directly  via  three  plane  wave  expansions 
according  to 


n(A2 -2/k2) 

2d~ d 
x  y 


t|/(i,i)  +  - 1:5 — j  iKi  „,i  h)  + - - — y 

2d  d  k2  a  b  2d  d  k2 
x  y  x  y 


by  invoking  PSIGEN  three  times. 

A  flow  diagram  for  ZMATG  is  given  in  Figure  15. 


3.3.2  PSIGEN 

Subroutine  PSIGEN  computes  the  i|>  defined  by  (49)  and  (36)  in  determin¬ 
ing  the  infinite  array  moment  matrices.  There  are  two  types  of  \|j  computations. 
One  type  sums  over  t|i  that  contain  both  positive  z  and  negative  z  traveling 

waves  as  defined  by  (36c).  The  other  sums  over  ^  that  contain  either  all 

P9 

positive  z  traveling  waves  or  all  negative  z  traveling  waves  but  not  both 
as  defined  by  (36a,b). 

Let  the  vector  w  „  be  the  difference  vector  between  one  of  the  end 
aS 

points  or  center  of  the  j  expansion  segment  and  one  of  the  end  points  or 
til 

center  of  the  i  test  segment.  The  five  possible  forms  for  w  are 
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Figure  15.  Flow  Diagram  for  ZMATG 
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Figure  15.  Continued 
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Figure  15.  Concluded 
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Similarly,  if 


|k(I  +  h./2)  ( * g_ > / 2  |  <  <  1 


(58) 


B+  is  approximated  by 


B.  = 


exp[-jk(wag  -  (i  ±  A  /2)f../2)-g_]  jk(I  +  A  12) 

—  -  ^  ^  J 


Equation  (57)  (or  (59))  can  be  derived  by  recasting  A+  of  (54)  (or  B+  of  (55)) 
into  the  form 

sinX 

X 

Under  the  condition  (56)  (or  (58)),  X  1  0  and  Y  r  A+  of  (57)  (or  B+  of  (59)). 

The  pq  looping  between  elliptical  contours  (circular  contours  in  recip¬ 
rocal  lattice  space)  is  described  in  the  description  for  subroutine  PSIA. 

A  flow  diagram  for  PSIGFF  is  given  in  Figure  16. 


3.3.3  PSIA 

Subroutine  PSIA  computes  Ha,  8)  of  (49)  with  the  if  (a,  6)  given  by 

P9 

(36  a, b)  .  The  plane  waves  of  (36a)  are  all  positive  z  traveling  waves,  those 

of  (36b)  are  all  negat  ive  z  travell ing  waves .  The  case  where  the  i(i  (a,  6)  are  de- 

P9 

fined  by  (36c)  (combination  of  positive  z  and  negative  z  traveling  plane 
waves)  is  handled  in  PS1GEN. 


In  PSIA,  the  if  (a,  8)  are  computed  from 
P9 


V  (~‘,i)  =  T7~. - tt — t— (exp[-jk(w  ,  -  £.)-g  ] 

PS  jkA.g  ( f  • g4 )  J  ap  2  j'  &±J 

J  2  J  - 


’xpl-jk(w^  +  -J-  r  ) -g4  J) 


(59) 
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ONE  PART  INTEGRATION 


Figure  16.  Flow  Diagram  for  PSIGEN 
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Figure  16,  Concluded 


for  wag.z  >  0  where  wflg 
If  |k(£  -g±)  A  /2 |  <  <  1, 


is  determined  from  one  of  the  five  forms  of  (52). 
however  iji  (a,  8)  is  computed  according  to  the  approx- 

ph 


imat ion 


(a,B)  -  exp[-jkw-g  +  ]  (60) 

PH  8Z 

The  doubly  infinite  pq  summations  in  (49)  are,  of  course,  approximated 

by  finite  summations.  The  area  in  pq  space  containing  the  points  over  which 

the  summation  is  performed  is  chosen  to  be  elliptical.  The  reason  for  chosing 

ellipses  is  that  a  characteristic  factor  of  the  li  functions  has  the  form 

pq 


exp[-jkR-g±] 


for  Rz  ?  0  where  R  is  a  spatial  vector  with  z  component  R^ .  For  p  and/or  q 
sufficiently  large,  g+  is  approximated  by 


g±  :  Px*  +  pyy  +  jP 


where 


x 


P 


y 


s  +  q 

y 


d 

y 


Then 

exp[-k | R  | p ] 

|F|  ;  - 1 - 

P 


Contours  of  constant  |f|  are,  therefore,  circles  of  radii  ■,  in  t.  ,  p  space 

x  y 

and  ellipses  in  p,q  space.  (The  space  defined  by  the  coordinates  and  r 

x  y 

is  typically  called  "inverse  lattice  space".) 
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The  summations  in  (49)  are  performed  by  adding  contributions  in  p,q 

space  corresponding  to  the  area  between  circular  contours  in  p  ,  p  space. 

x  y 

A  maximum  radius  p  is  chosen  bv  the  user,  and,  to  monitor  convergence,  a 
max 

"minimum"  radius  p  .  and  incremented  radius  A 
mm  p 


also  are  chosen.  The  sum¬ 


mation  over  the  is  performed  by  first  including  all  values  corresponding 

to  points  within  the  p  ^  circular  contour.  The  contribution  from  points 

within  the  p  .  and  p  .  +  Ap  circles  then  are  added  to  the  summation  and  so 

min  mm 

on  until  the  contribution  corresponding  to  all  points  within  the  p  circle 

max 

is  determined.  The  value  of  -p  at  each  step  in  the  process  can  be  recorded 
as  a  means  for  observing  convergence.  For  each  step  in  the  process,  the 
elliptical  contours  in  pq  space  that  correspond  to  the  circular  contours  in 


x'  1  y 

takes  on  the  values  p  =  p  .  ,  p  .  +  Ap ,  p  .  +  2Ap,  ,  p  and,  for  each 

min  ’  min  ’  mm  max 

p,  p  and  q  loops  where  the  bounds  on  p  and  q  are  the  ellipses  corresponding 


to  circles  in  p 


p^  space  of  previous  o  value  and  present  p  value  radii. 


Large  values  of  are  required  for  small  values  of  k  | ] 


When 


where  a. 

1 

is  used. 


t  h 

=  radius  of  j  segment,  the  user  supplied  value  of  p  (p 

max  max 

Otherwise,  p  is  computed  in  PSIA  according  to 
max 


Pu) 


p 


max 


(61) 


The  parameters  p  ^  and  Ap  also  are  adjusted  if  necessary,  and  Pmax  is  set 

to  a  minimum  value  (usually  =  5)  if  p  ,  as  computed  by  (61),  should  be  less 

max 

than  5. 

The  <!>  computation  in  PSIGEN  always  uses  p  =  Pu  since  the  character¬ 
istic  factor  is 


and  the  series  in  (49)  is  absolutely  divergent.  The  terms  in  the  series 
oscillate,  however,  and  the  series  appe>ars  to  converge  conditionally. 
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A  flow  diagram  of  PSIA  is  given  in  Figure  17. 


3.3.4  OFFD 


The  ij  off-diagonal  element  of  the  ground  screen  absent  moment  ma¬ 
trix  [Z]  or  the  ijt'1  element  of  the  ground  screen  "image"  matrix  [Z^ma®e]  is 
computed  in  OFFD.  The  equation  that  is  implemented  is  discussed  in  the  ZMATG 
description.  A  flow  diagram  of  OFFD  is  given  in  Figure  18. 


3.3.5  GNDS 

Subroutine  GNDS  computes  the  matrix  [Z^111386].  The  ijC^  element  of 
[Zima®e]  is  the  interaction  between  the  segment  ground  screen  image  and 
the  i*"^1  segment.  Toeplitz  symmetry,  even  applied  to  a  straight  wire,  is  not 
generally  exhibited  in  [Zlma8e],  Thus  all  ij  elements  are  computed  indivi¬ 
dually.  The  generally  wide  separation  between  wire  segments  and  images, 
however,  results  in  rapid  convergence  of  (49)  and  fast  computation  times.  A 
flow  diagram  is  given  in  Figure  19. 

3.3.6  VIDCJ 

This  subroutine  computes  the  generalized  voltage  excitation  column 
vector  V3a  (or  V^  )  of  the  ith  element  of  an  array  due  to  the  jth  neighbor¬ 
ing  radiating  element  current  distribution  in  the  region  of  influence.  Each 
segment  of  the  j*'*'  element  current  distribution  is  approximated  as  a  Hertzian 

—fla  -bb 

dipole.  The  voltage  is  accumulated  in  (or  )  as  VIDCJ  is  called  for 
each  radiating  element  in  the  region  of  influence. 

A  flow  diagram  of  VIDCJ  is  given  in  Figure  20. 


3.3.7  AOFI 

This  subroutine  computes  region  of  influence  boundaries  for  the  i*'*' 
radiating  element  in  a  higher  approximation  computation.  All  radiating  ele¬ 
ments  within  the  region  of  influence  will  be  included  in  the  higher  approxi¬ 
mation  computation  of  Vaa  (or  V^^)  corresponding  to  the  it^'  element.  Lens 
face  elements  are  tested  for  location  within  the  region  of  influence  by  loop¬ 
ing  over  all  cells  and  considering  each  cell  array  description  individually. 
All  array  descriptions  are  considered  "regular"  for  this  computation, 
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BEGIN 


Figure  17.  Flow  Diagram  for  PSIA 


66 


END 

END  LOOP 


Figure  19.  Flow  Diagram  for  GNDS 


Figure  20.  Flow  Diagraa  for  VIDCJ 
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i.e.,  the  array  boundaries  are  parallel  to  the  cell  lattice  coordinates. 
Judicious  choice  of  IARY  (L,2)  by  the  user  when  defining  cells  by  array  de¬ 
scriptions  (as  done  for  FAM)  is  called  for.  Note  that  for  irregular  arrays 
IARY  (L,3)  4  IARY  (L,5),  the  user  may  choose  IARY  (L,2)  equal  to  neither 
IAP.Y  (L,3)  or  IARY  (L, 5).  Since  only  IARY  (L,2)  is  used  by  AOFI  in  deter¬ 
mining  the  region  of  influence,  several  elements  in  an  irregular  array  would 
not  be  tested  for  location  within  the  region  of  influence  (an  irregular  ar¬ 
ray  is  "regularized").  The  user  must  be  aware  of  this. 


A  flow  diagram  of  AOFI  is  given  in  Figure  21. 


3.3.8  CELPAT 

I 

This  subroutine  computes- the  radiation  pattern  of  a  cell  described  in 
terms  of  several' arrays.  The  coding  follows  the  theory  of  [5,  Section  3.4] 
where  "regular"  and  "irregular"  shaped  arrays  are  discussed.  The  6  and  <f> 
(global  system)  components  of  J-field  that  results  are  normalized  by 
exp  [-jkR]/(4nR) .  ‘ 

A  flow  diagram  of  CELPAT  is  given  in  Figure  22. 

3.3.9  MUTUAL 

s.b 

The  array  b^  to  array  a  ground  screen  penetration  mutual  impedance  z 
is  computed  in  subroutine  MUTUAL.  The  induced  short  circuit  current  distri¬ 
bution  on  the  array  ja  reference  element  aiso  is  computed.  The  coding  follows 
closely  the  theory  in  the  latter ’part  of  Section  2.2.  The  good  conducting 
screen  approximation,  -  -,z ,  is  employed.  A  flow  diagram  is  given  in 
Figure  23. 

3. 3. 10  ZIMPS  * 

Subroutine  ZIMPS  computes  the  moment  matrix  for  an  isolated  thin 
wire  radiating  element  over  a  perfect  ground  screen.  This  matrix  is  needed 
in  the  higher  approximation  method  (HAM).  An  option  for  removing  the  ground 
screen  from  the  computation  exists  by  setting  IGSC  “  0.  The  user  must  not 
exercise  this  option,  however,  in  a  HAM  computation. 

The  thin  wire  moment  method  code,  upon  which  ZIMPS  is  based,  is  the 
code  WRSMOM  available  from  NAPS  [7].  Since  WRSMOM  assumes  geometry  data  in 


i 
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Figure  21.  Flow  Diagram  for  AOFI 


COMPUTE  DOT  PRODUCT  BETWEEN  CURRENT  SEGMENT  POSITION 
VECTOR  AND  FIELD  POINT  VECTOR  AND  EXPONENTIAL  THEREOF 


:o: 


Figure  22.  Continued 
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Figure  22.  Concluded 


BEGIN 


Figure  23.  Flow  Diagram  For  MUTUAL 
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m 


Figure  23.  Concluded 
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wavelengths,  this  data  is  converted  to  wavelengths  in  ZIMPS.  The  matrix  ele¬ 
ments  then  are  computed  by  calls  to  subroutines  ZIMP2  or  ZIMP3,  A  flow  dia¬ 
gram  is  given  in  Figure  24 

3.3.11  LENSIN,  CELLIN.  ARAYIN,  GEOMET,  TPIN 

These  subroutines  control  the  input  data  defining  the  lens  structural 
and  electrical  characteristics,  lens  illumination,  software  control  param¬ 
eters,  and  radiation  pattern  plot  format.  A  user  interactive  manner  is  es¬ 
tablished  whereby  the  program  queries  the  user  for  each  data  input. 


Subroutine  LENSIN  requests: 

•  existence  of  a  radiation  pattern  file.  If  "yes",  the  logical 
variable  FILFG  is  set  equal  to  TRUE,  and  pattern  data  computed 
is  superimposed  on  the  existing  pattern  data. 

•  direction  cosine  limits  of  pattern  calculation.  If  "12",  the 

y  axis  direction  cosine  will  vary  most  rapidly.  If  "21",  the 

x  axis  direction  cosine  will  vary  most  rapidly.  Input  of  the 

x  axis  direction  cosine  minimum,  maximum,  and  increment  are 
followed  by  those  for  the  y  axis  direction  cosine. 

•  boom  length. 

•  ground  screen  transmission  coefficients  for  E-field  perpendicu¬ 
lar  and  parallel  to  plane  of  incidence. 


Subroutine  CELLIN,  for  each  cell,  calls  ARAYIN  and  GEOMET  and  also 
requests: 


•  cell  ID  number,  number  of  plane  waves  in  illumination,  and  num¬ 
ber  of  feedback  inodes  to  be  considered  in  analysis. 


reciprocal  lattice  space 
incremental  radii  (p  , 


(p  ,  p  space)  minimum,  maximum,  and 
x  y 

p  ,  Ap)  used  in  the  i|/(a,B)  summations, 
max 


•  cell  coordinate  system  origin. 

•  cell  coordinage  system  rotation  relative  to  global  system. 

•  lattice  parameters  d  ,  d  ,  A  . 

x  y  y 

•  cell  area  via  array  description  (call  to  ARAYIN). 

•  characteristics  of  modules  and  feedlines  to  form  limited  scat¬ 
tering  parameter  or  impedance  parameter  models. 

•  array  a.  and  array  b  radiating  element  data  (call  to  GEOMET). 


Subroutine  ARAYIN  requests: 


•  number  of  array  descriptions. 

•  each  array  description  in  terms  of  number  of  columns  (elements 
along  line  parallel  to  lattice  e  coordinate),  number  of  rows 
(elements  along  line  parallel  to  y  coordinate),  number  of  ele¬ 
ments  in  first  column  (set=0  if  array  is  regular  [5,  Section  3.4]) 


77 


Figure  24.  Flow  Diagram  for  ZIMPS 
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Figure  24.  Concluded 


number  of  elements  in  last  column,  and  sum-subtract  code  =  1 
or  -1  (not  used  in  present  versions  of  FAM  or  HAM) . 

•  cell  coordinates  of  array  reference  element  (must  be  the  "lower 
left  hand  corner"  [5,  Figure  3-16]). 

•  array  lattice  parameters  (may  differ  from  cell  lattice  param¬ 
eters  d^,  d  ,  A  when  describing  elements  for  which  a  higher 

approximation  analysis  is  sought  (HAM)). 

Subroutine  GEOMET  requests: 

•  frequency. 

•  folded  dipole  model  control  parameter  ("yes"  or  "no"1'.  If 
"yes",  then  folded  dipole  length  and  characteristic  impedance 
of  transmission  line  mode  are  entered. 

0  radiator  geometry  as  a  collection  of  straight  wires  (fo'Jov 
closely  [71).  The  user  provides  the  number  of  wires  and.,  ‘‘or 
each  wire,  beginning  and  end  point  coordinates,  radius,  r ■.•••• 
nection  code,  number  of  feed  points  (0  or  1  (  present  vi¬ 
sion  of  FAM  -  one  feed  point  pet  radiator),  and  location  of 
feed  points.  The  user  must  take  care  that  the  numbers  of 
segments  per  wire  and  wire  lengths  yield  nearly  uniform  seg¬ 
ment  lengths  throughout  any  noiuii&Ji int  collection  of  wires. 
Segment  length  variations  greater  than  ^1.5:1  are  discouraged. 
Also,  the  algorithm  the  program  uses  to  segment  a  wire  depends 
somewhat  on  connection  code.  It  1,  is  the  wire  length,  N  the 
number  of  segments,  and  A  the  segment  length  then: 

no  connection:  A  -  L / ( N  +1) 

s 

one  connection:  A  =  l./(N^  +  0.5) 

two  connections:  A  =  l./N 

s 

The  connection  code  for  a  wire  that  touches  the  ground  screen 

must  indicate  a  connection  at  that  end  of  the  wire. 

Subroutine  TPTN  request,  for  T-plane  radiation  pattern  plot: 

0  x  axis  direction  cosine  limits. 

0  y  axis  direction  cosine  limits. 

0  minimum  level  of  ffeld  value  to  be  assigned  a  distinguishing 
plot  notation  ("FLOOR")  and  field  level  increment  between 
distinguishing  notations. 

Subroutine  CF.LI.1N  computes  module  progressive  phase  settings  to  ac¬ 
commodate  beam  scanning  and  focusing.  Subroutine  GEOMET  computes 
wavelength,  wave  number,  reference  element  segment  position  vectors 
(to  the  segment  centers),  and  segment  unit  vectors. 


3.3.12  FEED 


Subroutine  FEED  computes  the  plane  wave  illumination  of  a  cell 
according  to 


-  -Jkr  n  P  ^2 

E  =  -  (-—) 

r  Att 


(P  x  +  P  y  +  P  z) 
x  y  z 


where 


n 

P 

e 

L 


y 


are  the  global  coordinates  of  the  cell  coordinate  origin 
(date  input) 

free  space  wave  impedance 
effective  radiated  power  (data  input) 


*  boom  length  (data  input) 


r  *  /x2+y2+L2 
g  g 

P  ,  P  =  x  and  y  polarization  components  (data  input) 
x  y 

■  z  polarization  component  computed  according  to 


P 

z 


P  s 

X  X 


+  P  s 
_ LJL 


to  assure  polarization  orthogonal  to  propagation  direction.  Here, 

s  are  the  propagation  direction  cosines  computed  according  to  s  ■ 
z  x 


s  ,  s  , 
x  y 

Vr’ 


s  ■  y  /r,  s  -  L/r. 

y  g  z 


3.3.13  SCAN 

Subroutine  SCAN  computes  the  scan  (active)  impedance  of  an  array 

over  a  perfect  ground  screen.  Here,  the  array  moment  matrix  [Z®8]  is  Invert- 

ps  —  X 

ed  via  LINEQ,  the  element  y  *  [Z"  corresponding  to  the  feed  port  seg- 

S’  S  -1 

ment  ig  is  extracted,  and  the  active  impedance  is  g4-'en  by  x  ■  y 


3.3.14  SCIA,  EX 

Subroutine  SCIA  computes,  for  a  cell,  the  array  a  reference  element 
short  circuit  current  distribution.  All  array  a  ports  are  short  circuited 
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and  the  ground  screen  is  present.  The  excitation  is  the  plane  wave  illumi¬ 
nation  as  defined  in  FEED. 

Subroutine  EX  is  called  by  SCIA  to  compute  the  voltage  excitation, 

column  vector  V  ,E>  according  to  (40).  The  reference  element  short  circuit 

current  distribution  is  then  computed  by  solving  (41)  with  V  removed.  The 

matrix  [Z®s]  has  been  inverted  in  SCAN  via  LINEQ.  The  reference  element 

short  circuit  port  current  1^  is  extracted. 

s 

3.3.15  PERP 

th 

Subroutine  PERP  computes  the  i  segment  radius  vector  a^  according 
to  the  discussion  following  (38).  It  the  segment  is  oriented  normal  to  the 
array  plane  ( £ ■  ±z)  the  double  cross  product  relation  for  is  not  defined. 

The  a^  then  is  chosen  to  be  a^  *  a^x. 

3.3.16  ZCOL 

Subroutine  ZCOL  extracts  a  specified  column  from  a  square  matrix. 

3.3.17  LINEQ 

Subroutine  LINEQ  inverts  a  complex  matrix. 

3.3.18  CEXP ,  TAN 

These  function  subprograms  compute  eX  where  x  is  complex  (CEXP)  and 

tan  (0)  where  6  is  in  radians  (TAN).  If  the  real  part  of  x  i  -70,  eX  is  set 

32 

to  zero.  If  cos  (f»)  «  0,  tan  (0)  -is  set  to  10 

3.3.19  ZIMP2 ,  2IMP3 ,  HAMPSI 

Subroutines  ZIMP2  and  ZIMP3  and  function  HAMPSI  are  minor  modifica¬ 
tions  of  coding  in  the  isolated  element  moment  method  program  WRSMOM  (7). 

3.3.20  XYOCS 

Subroutine  XYCCS  transforms  x,  y  cell  coordinates  to  lens  system 
(global)  coordinates. 

3.3.21  TPLANE 

Subroutine  TPLANE  plots  the  radiation  pattern  versus  x  and  y  direc¬ 
tional  cosines.  Printer  characters  are  used  to  denote  field  levels.  This 
routine  was  written  by  R.  Hancock  [8]. 
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SECTION  4 
RESULTS 


The  pulse  expansion-point  matching-finite  difference  operator  moment 
method  derived  in  Section  2.2  was  applied  to  the  infinite  rectangular  lattice 
array  of  straight  halfwave  dipoles  previously  analyzed  by  Chang  [9).  The 
lattice  parameters  were  d^  =  0.25X,  d^  =  1.2X,  and  Ay  =  0  (rectangular  lat¬ 
tice).  The  E-plane  and  H-plane  scan  active  admittances  are  shown  in  Fig¬ 
ure  25.  The  excellent  agreement  compares  well  with  the  previous  comparison 
between  Chang's  results  and  the  pulse  expansion-pulse  matching  moment  method 
originally  employed  in  the  lens  simulator  [5,  Figure  5-3].  The  finite  dif¬ 
ference  operator  method  is  considerably  more  efficient  also:  only  13  seg¬ 
ments  subdivided  a  dipole  for  the  finite  difference  results  whereas  21  were 
required  for  the  pulse  matching  results.  A  Pmax  of  100  (Section  3.3.3)  were 
used  in  obtaining  the  Figure  25  results. 

Active  Impedances  and  current  distributions  of  infinite  arrays  of 
straight  and  "swept  back”  dipoles  on  equilateral  triangular  lattices  also 
were  computed  using  the  lens  simulator  finite  difference  moment  method  dis¬ 
cussed  in  this  report.  Figure  26  shows  a  straight  halfwave  dipole  array  ra¬ 
diating  element  (y  polarized)  a  quarter  wavelength  in  front  of  a  perfect 
ground  screen.  All  dipoles  were  assumed  to  be  fed  through  perfect  baluns. 

The  radiating  nature  of  the  feedlines,  then,  was  that  of  parasitic  scatterers 
not  physically  connected  to  the  dipole  feed  ports.  The  two  types  of  feedline 
scatterers  modeled  are  shown  in  Figure  26.  One  was  a  0.2X  straight  feedline 
and  the  other  a  0.275X  bent  feedline.  Both  feedlines  lay  parallel  to  the 
yz  plane  (coordinate  system  of  Figure  2  applies).  All  array  elements  were 
located  on  the  equilateral  triangular  lattice  defined  by  Figure  4  with 

d  -  0.715X 
x 

d  -  0.826X 

y 

Ay  -  0.413X 

The  H-plane  scan  active  impedance  for  the  array  without  feedline  scatterers 
is  plotted  in  Figure  27.  The  effect  of  the  H-plane  grating  lobe  singularity 
occurring  at  23.49°  is  evident.  The  ground  screen  removes  the  infinity  that 
would  otherwise  occur  at  an  H-plane  grating  lobe  singularity  for  straight 


T  - *■-- 
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E-PLANE  REAL 


SCAN  ANCLE  (DEC) 


ri«ure  25.  Reference  Element  Active  Admittance 
for  a  X / 2  Dipole  Array  0.25A  Above 
a  Perfect  Ground  Screen.  The  Dipole 
Radii  Are  0.007022A.  The  Lattice 
Parameters  Are  d  ■  0.25X,  d  =  1.2A, 
Ay  -  0  *  y 
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dipoles.  A  similar  plot  for  straight-dipole  array  with  straight  feedline 

scatterers  differed  very  little  from  Figure  27.  The  absence  of  a  stronger 

singularity  effect  at  23.49°  scan  likely  is  a  consequence  of  H-plane  scanned 
■+ 

E-field  radiated  by  the  dipoles  along  the  array  plane  being  cross  polarized 
with  respect  to  the  feedlines. 

The  E-plane  scanned  active  impedances  for  the  straight-dipole  array 
for  the  three  cases:  feedlines  absent,  straight  feedlines,  and  bent  feedlines, 
are  shown  in  Figure  28.  The  feedline  absent  and  bent  feedline  cases  appear 
strongly  effected  by  the  E-plane  grating  lobe  singularity  occurring  at  29.22°. 
The  straight  feedline  case  exhibits  rapid  impedance  variation  in  the  vicinity 
of  22°  scan.  The  reference  element  current  distributions  plotted  in  Fig' 
ures  29,  30,  and  31  also  bear  out  this  behavior.  The  reactive  current  on  the 
straight  feedline  swings  widely  from  capacitive  to  inductive  in  scanning 
from  20°  to  25°  (Figure  30).  The  bent  feedline  current  does  not  demonstrate 
this  behavior. 

E-plane  scan  results  for  the  Figure  32  dipole  array  without  feedlines 
and  with  straight  feedlines  are  shown  in  Figures  33  and  34.  Without  feedlines, 
the  swept  back  dipole  array  active  impedance  varies  less  with  F-plane  scan 
than  does  that  of  the  straight-dipole  array.  The  effect  of  straight  feed¬ 
lines  on  the  E-plane  active  impedance  is  at  least  as  pronounced  for  the 
swept  back-dipole  array  as  it  is  for  the  straight-dipole  array. 

A  pmx  of  50  (Section  3.3.3)  was  used  in  obtaining  the  equilateral 
lattice  results.  Sampled  computations  with  Pmax  =  100  were  within  ^10  per¬ 
cent  of  the  p  =50  results.  Also  9  segments  subdivided  the  dipoles  and 
max  ° 

5  segments  the  feedlines. 
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Figure  28.  E-Plane  Active  Impedances  of  Straight 
Dipole  Array  Without  Feedlines,  With 
Straight  Feedlines  and  With  Bent  Feed 
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Figure  29.  Straight  Dipole  Array  Reference  Element 
Current  Distribution:  1  V  Source,  No 
Feedlines 
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Figure  31.  Straight  Dipole  Array  Reference  Elejnent  Current 
Distribution:  1  V  Source,  Bent  Feedlines 
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STRAIGHT  FEEDLINE 


Figure  32.  Swept  Back  Dipole  Array  Element  Above  a 
Perfect  Ground  Screen  With  a  Straight 
Feedline  Scatterer 
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SWEPT  BACK  DIPOLE  (1) 


Figure  34.  Swept  Back  Dipole  Array  Reference 
Element  Current  Distribution:  1  V 
Source,  Mo  Feedline 
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